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With the experimental data from STAR, PHENIX and BRAHMS on the centrality and rapidity 
dependence of the pr spectrum in Au+Au and d+Au collisions, we show that there exists a scaling 
distribution which is independent of the colliding system, centrality and rapidity. The parameter for 
the average transverse momentum (pr) increases from peripheral to central d+Au collisions. This 
increase accounts for the enhancement of particle production in those collisions. A non-extensive 
entropy is used to derive the scaling function. 
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The distributions of produced particles play very im- 
portant role in the study of new state of matter produced 
in relativistic nuclear collisions. From the relationship 
between the distributions at different centralities and ra- 
pidities one can see whether there is a suppression or 
enhancement, thus dynamical mechanism for those par- 
ticle production can be revealed. In earlier works a 
scaling distribution for produced pions in Au+Au colli- 
sions in mid-rapidity region at RHIC energies was found 
that is independent of the centrality and collidingenergy. 
Similar scaling behaviors have been found in [2, 3]. The 
existence of such scaling behaviors of the spectrum for a 
species of particle may be an indication of universal un- 
derlying particle production mechanism in the collisions. 
One may ask two questions now. (1) Can the scaling 
be extended to non-central rapidity region? (2) Can a 
scaling behavior similar to that shown in [l[ be found for 
d+Au collisions at different centralities and rapidities? If 
the answers to above two questions are yes, one should 
investigate the difference between the scaling functions 
for different colliding systems. In this paper, we look for 
the scaling behavior of particle distributions in Au+Au 
collisions in noncentral rapidity region, in d+Au colli- 
sions at different centralities and rapidities at RHIC en- 
ergies. The same universal scaling behavior is found for 
both Au+Au and d+Au collisions at different centralities 
and rapidities. A possible origin for the universal scal- 
ing distribution is given from the principle of maximum 
non-extensive entropy. 



Generally, the invariant particle distribution E^^- for 
a given species of particles depends on the collision cen- 
trality (3, (pseudo)rapidity y (or rf) and the transverse 
momentum px- When we focus on the particle produc- 
tion in a collision with given centrality in a small rapidity 
region around y, we should distinguish two different fac- 
tors. One is the number density dN/dy in rapidity, which 
tells us how many particles are produced in the region. 
The other is the probability for a produced particle hav- 
ing transverse momentum px- It is this probability that 
dictates the shape of the transverse momentum distri- 
butions. In this paper, we are interested only in such a 
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FIG. 1: Scaling distribution 9(u) showing the coalescence of 
four rapidity bins in the BRAHMS data for h~ production 
at 7j = 2.6 and 3.2 in Au+Au collisions. The data points are 
taken from Fig.l in Q. 



probability which can be defined as 



<P( PT ,y,f3) = 2TT / — , 



,d 3 N I dN 
dy 
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where 



<1N 



27T J E^fpxdpx is the multiplicity of par- 
ticles produced in a unit rapidity region. So Q(px,y,(3) 
satisfies the normalization condition 



J ®{PT,y,P)PTdp T = f • 



(2) 



In the definition of the centrality and rapidity de- 
pendence is written explicitly. Experimental data show 
similar shape of the distributions for different centrali- 
ties and rapidities. Because of the similarity in the shape, 
those distributions of pions produced in Au+Au collisions 
in the mid-rapidity region at difference centralities have 
been put together in [l|] to the same curve by shrinking px 
by a factor K and shifting the distribution by another fac- 
tor A. Then a scaling variable is defined as u = px / (pt), 
and the normalized distribution is obtained as a function 
^(u) = (px) 2 $((px}u,y, (3). Here the average transverse 
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momentum is defined as 

(p T ) = J <S>{p T ,y,P)PTdpT ■ (3) 

The scaling function reads 

= 2.1 x 10 4 (u 2 + 11.65)~ 48 [1 + 25exp(-1.864u)].(4) 

The obtained ^(u) is universal in the sense that it is in- 
dependent of the centrality (3 at mid-rapidity in Au+Au 
collisions. 

Now we extend the discussion to the particle produc- 
tion in non-central rapidity region in Au+Au collisions 
at 200 GeV. For this purpose, we work on the experimen- 
tal data Q for negatively charged particles hr at large 77 
with exactly the same procedures as in [l[ and try to put 
all the data points to the curve for the scaling function 
fy(u). It turns out that the scaling function given by Eq. 
(j4]) can also describe the charged particle production at 
77 =2.6 and 3.2 for different centralities, as can be seen in 
Fig. [TJ The shrinking factor for pt is given in TABLE fl] 
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TABLE I: Fitted parameter (pr) in the scaling distribution 
for Au+Au and d+Au collisions at ^Jsnn = 200 GeV, for 
different centralities and rapidities shown in Figs. (1-3). 

The second observation in this paper is that the same 
scaling pr distribution can describe the pr distributions 
of charged particles produced in d+Au collisions [5(] and 
7T° data Qin central rapidity region at different centrali- 
ties, as shown in Fig. [2] The fitted values of parameter 
(p T ) for the PHENIX data are also tabulated in TABLE 
I. 

To investigate the rapidity dependence of the distribu- 
tion in d+Au collisions, we take the data for n + and ir~ 
at y — 3 for different centrality cuts 0] and try to put 
the data points to the obtained scaling curve with exactly 
the same procedure as mentioned above. The agreement 
is quite good, as shown in Fig. [31 The fitted values of 
the parameter (pt) are also given in TABLE I. Then we 
can conclude that the same scaling function can describe 
pion and charged particle production in both Au+Au 
and d+Au collisions at different centralities and rapidi- 
ties. In this sense, we can say that the obtained scaling 
distribution is universal. 

To show the degree of agreement between the exper- 
imental data and the fitted results, we define B as the 




FIG. 2: Scaling distribution ty(u) showing the coalescence 
spectra from the STAR data for charged particle and from 
PHENIX 7T° data production in d+Au collisions at mid- 
rapidity. The data points are taken from [1,0]. 
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FIG. 3: Scaling distribution ^(u) showing the coalescence of 
four centrality bins in the BRAHMS data for tv~ production 
in d+Au collisions at forward direction. The data points are 
taken from Q. 

ratio of the data to the fitted results. For illustration, 
B, calculated from results in Fig. 3, as a function of u 
is shown in Fig. 4. Obviously, values of B are very close 
to 1, indicating that the scaling of the produced pion 
spectrum is true with high accuracy. 

Once we have found a scaling law for the distributions 
for Au+Au and d+Au collisions for different centrali- 
ties and rapidities, we need to dig out the physical im- 
plications behind the scaling behavior. From the scal- 
ing function ^(u), one can see that the only parame- 
ter characterizing the normalized distribution is the av- 
erage transverse momentum (j>t) which depends on the 
centrality, rapidity and the colliding energy and system. 
Once (pt) is known, both the soft part with low pt and 
the hard part with high pt are determined from ^(u). 
Traditionally it was assumed that the bulk of particles 
is produced by soft processes which determine (pt)- If 
the soft and hard processes can be separated, as used in 
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FIG. 4: Ratio B of the experimental data and the fitted re- 
sults, as shown in Fig. 3, for the pion production as forward 
direction. 



the two-component model [8(, and each part has differ- 
ent dynamical origin, it is extremely hard to imagine that 
they can be united into a single scaling function. Thus, 
we conclude that the hard process must contribute con- 
siderably to the average px and is a coherent part of 
particle production. In fact, we have known from the 
study of particle correlations in central Au+Au collisions 
at RHIC energies that the lost energy from hard partons 
must heat the medium and increase the mean transverse 
momentum. 

On the other hand, we have known from experiments 
that the particle production is suppressed, called jet 
quenching 0] , in central Au+Au collisions because of the 
possible formation of dense hot medium but enhanced in 
d+Au collisions due to the Cronin effect [13], possibly 
because of the px broadening from multiple soft colli- 
sions before a hard one or from the recombination 
of soft partons with those from showers [12]. One can 
show that the enhanced production in d+Au collisions 
can be accounted for from our fitted results that the mean 
transverse momentum (px) increases from peripheral to 
central collisions. We can define a new variable Qcp to 
measure the suppression or enhancement from the scaling 
function $(pr,y, 0) as 



Qcp = 



<§>(p T ,y,P)c _ (^(pt/(pt))/(pi 
®(PT,y,P)p 



(*(pt/(pt))/(pt) 2 ) j 



(5) 



An advantage of this definition of Qcp over the tradi- 
tional Rcp is that it involves only the probability for a 
produced particle having certain transverse momentum 
and is independent of the total multiplicity. So the over- 
all enhancement in the backward region (y < 0) in d+Au 
collisions cannot be shown in Qcp- bi calculating the 
value of Qcp we do not need to use the number of binary 
hard collisions, which is model dependent. Qcp is shown 
in Fig. [5] as a function of px for d+Au collisions at mid- 
rapidity and it increases slowly to about 1.33. Thus the 



production of high px pions is enhancement by about 30 
percent. This is in agreement with experimental obser- 
vations [13] ] . In this way the enhancement or suppression 
behavior can be shown clearly in a model independent 
way. 
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FIG. 5: Qcp as a function of pr for d+Au collisions, 
centrality cuts are 0-20% and 40-100%, respectively. 
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Finally, we investigate nucleon-nucleon collisions at 
■s/s = 200 GeV to sec whether the particle spectrum can 
be described by the same scaling function. The experi- 
mental data on the pr distribution in the central rapidity 
region in p+p collisions can be parameterized [Tij as 



dN 
p T dp 7 



A(a+p T y 



(6) 



where A, a and n can be determined from fitting the data. 
It is astonishing that this parameterization for pion dis- 
tribution in p+p collisions, with the parameters suitably 
chosen, can also be put to the scaling curve for ^(it) by 
shrinking horizontally and shifting vertically. As shown 
in Fig. 6 the n + spectrum in p+p collisions can be fit- 
ted to the scaling function. Then we come to our final 
conclusion that the same scaling function can be used, 
for different colliding systems at different centralities and 
rapidities, to give the probability for a produced particle 
having px once the parameter (px) for the average trans- 
verse momentum is known, thus the scaling function is 
universal. 

Since ^(u) is universal, the centrality and rapidity 
dependence in &(px,y,(3) is totally encoded in (px)- 
Straightforwardly, we can have 



(7) 



Therefore, the ratio (p T ) / ' (px) n is a constant independent 
of the colliding system, the centrality, and rapidity. This 
consequence can be easily checked experimentally. For 



(Pt) = J PT$(PT,y,f3)PTdp T 
= (p r )" / u n ^(u)udu . 
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n = 2, 3 and 4, the ratio equals to 1.65, 4.08 and 14.39, 
respectively. 
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Eq. (4) 




----Eq.(9) 
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FIG. 6: Comparison between from fitting the data and 

from the maximum entropy method. The insert gives the 
ratio between them. The points are ir + distribution in p+p 
collisions taken from [l5l ]. 

We are now in a position to derive an expression for the 
particle distribution from some statistical consideration. 
One method is from the principle of maximum entropy. 
The scaling function is not an exponential form, therefore 
the well-known Shannon entropy cannot be used for our 
purpose. Instead, as suggested in [hq , a non-extensive 
entropy may be used, considering the complex nonlinear 
interactions and the possible presence of fractal chaotic 
behavior in nuclear collisions, to derive an expression for 
the particle distribution. The entropy we will use here is 
the Tsallis entropy [17] which reads for our case 



S q = (l- / &(p T ) PT dp T )/(q-l) , 



(8) 



with parameter q as a measure of the non-extensiveness of 
the system. The maximum entropy under conditions de- 
manded by Eqs. and ([3]) gives, after the same shrink- 
ing and shifting as above, 



= A(a + u)~ 



(9) 



withn = l/(l-q) and a = (n-3)/2,A= a n ~ 2 (n-l)(n- 
2). So there is one parameter q to determine the shape 
of the distribution. The functional form of 'J'(w) derived 
from the maximum entropy method is very different from 
that given by Eq. J4} but is the same as that used in pa- 
rameterizing the p+p data. Then one can compare the 
distributions ^(u) from our fitting procedure and from 
the maximum entropy. As shown in Fig. [SI they agree 
with each other extremely well in a very wide range of 
u when n — 13.2. There in log scale one can hardly see 
any difference between them. Small difference can be 
seen if we show a ratio between them in a linear scale. 
As can be seen from the insert in Fig. 6 the largest dif- 
ference is less than ten percent, though the distribution 



covers six orders of magnitudes in the shown range of 
u. This agreement indicates that the produced particle 
system, including both the soft and hard particles, may 
be a non-extensive fractal system and the fractal prop- 
erty is independent the colliding system. The parameter 
characterizing the fractal property of the system is then 
q = 1 — 1/n = 0.924, not too far from 1. 

In summary, a scaling distribution for pions is found 
for Au+Au, d+Au and p+p collisions at RHIC energy. 
The scaling distribution is universal, since it does not 
depend on the colliding system, the centrality and ra- 
pidity. What characterizes the dynamical difference in 
the particle production is the value of the average trans- 
verse momentum (pr). A simple prediction which can 
be checked experimentally is that (Pt)/{Pt)" is indepen- 
dent of the colliding system, the colliding centrality and 
rapidity. Possible statistical origin for the scaling distri- 
bution of pions is discussed from a non-extensive entropy. 
Different particles may have different scaling properties. 

This work was supported in part by the National 
Natural Science Foundation of China under grant No. 
10475032. 



[2 
[3 

[4; 

[6 

[-: 

[8 

[9 
[10 

[11 

[12 

[13 

[14 
[15 

[16 
[17 



Rom. Rep. Phys. 
Lett. 86, 3496 



R.C. Hwa and C.B. Yang, Phys. Rev. Lett. 90, 212301 

(2003) ; Phys. Rev. C 67, 034902 (2003); Phys. Rev. C 
67, 064902 (2003). 

I. Zborovsky and M. Tokarev, hep-ph/060322l| and ref- 
erences therein. 

J. Schaffner-Bielich, D. Kharzeev, L. McLerran, and R. 
Venugopalan, Nucl. Phys. A 705, 494 (2002). 
C. Ristea (for BRAHMS Collaboration), AIP. Conf. 
Proc. 828, (2006) 353. 

J. Adams et al., (STAR Collaboration) , |nucl-ex704 08016 
Phys. Rev. C 70, 064907 (2004). 

S.S. Adler et al., PHENIX Collaboration, 
|nucl-th/0610036| 

H. Yang (for BRAHMS Collaboration) 
58, (2006) 5. 

X.N. Wang and M. Gyulassy, Phys. Rev 
(2001). 

X.N. Wang, Phys. Rev. C 58, 2321 (1998). 
J.W. Cronin et al., Phys. Rev. D 11, 3105 (1975); D. 
Antreasyan et al, Phys. Rev. D 19, 764 (1979 ). 
For a review, see A. Accardi, hep-ph/0 2 1 2 1 48 1 and refer- 
ences therein. 

R.C. Hwa and C.B. Yang, Phys. Rev. Lett. 93, 082302 

(2004) . 

S.S. Adler et al., PHENIX Collaboration, Phys. Rev. C 
74, 024904 (2006). 

R. Hagedorn, Riv. Nuovo Cimento 6, 1 (1983). 

J. Adams et al., (STAR Collaboration), |nucl-th/0601033| 

Phys. Lett. B 637, 161 (2006). 

A.S. Parvan, T.S. Biro, Phys. Lett. A 340, 375 (2005). 
C. Tsallis, J. Stat. Phys. 52, 479 (1988); Braz. J. Phys. 
29, 1 (1999). 



